Gapped Two-Body Hamiltonian for continuous-variable quantum computation 
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We introduce a family of Hamiltonian systems for measurement-based quantum computation with 
continuous variables. The Hamiltonians (i) are quadratic, and therefore two body, [ii) are of short 
range, {Hi) are frustration-free, and (iv) possess a constant energy gap proportional to the squared 
inverse of the squeezing. Their ground states are the celebrated Gaussian graph states, which are 
universal resources for quantum computation in the limit of infinite squeezing. These Hamiltonians 
constitute the basic ingredient for the adiabatic preparation of graph states and thus open new 
venues for the physical realization of continuous- variable quantum computing beyond the standard 
optical approaches. We characterize the correlations in these systems at thermal equilibrium. In 
particular, we prove that the correlations across any multipartition are contained exactly in its 
boundary, automatically yielding a correlation area law. 
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Introduction. — The realization of a device that can per- 
form arbitrary quantum-state manipulations - a univer- 
sal quantum processor - is nowadays one of the most 
active and promising searches in physics. All the re- 
sources on vifhich we count for the realization of these 
machines in turn fall into one out of two fully-general cat- 
egories: quantum states of discrete- variable (DV), finite- 
dimensional systems and quantum states of continuous- 
variable (CV), infinite-dimensional ones. For both cate- 
gories, considerable progress has been achieved. In par- 
ticular, an entire spectrum of new possibilities for state 
manipulation was opened by the landmark discovery that 
it is possible to process quantum information in a uni- 
versal manner by the simple act of performing local 
measurements [1]. This discovery originally took place 
in the finite-dimensional scenario, but was later on ex- 
tended to CV systems [5]. In these measurement-based 
quantum-computation (MBQC) approaches, information 
processing proceeds by a sequence of adaptive single- 
particle measurements on massively entangled multipar- 
ticle states prepared in advance. These states are the 
so-called cluster states, introduced first for DV systems 
[HIS] and then extended to the CV case Hg]. The local 
measurements consume cluster-state entanglement as the 
main resource of the computation. 

Cluster states are in turn particular instances of a more 
general family: the graph states [5]. Other examples of 
graph states of importance are the Greenberger-Horne- 
Zeilinger states and many codewords for quantum error- 
correction [6] . Every graph state is associated to a math- 
ematical graph G'(v,£) = {^i^}, of vertices i G V and 
edges {«, j} e £, with 1 < i, j < N. In both the DV and 
CV cases, graph states can be created in an operational 
way: starting from a product state of N particles, each 
one associated to a vertex in V, a sequence of entangling 
operations between every pair of particles connected by 
an edge in £ is applied to obtain the desired state [5] . 



Alternatively, for DV systems, there exists a conceptu- 
ally different approach: Every finite-dimensional graph 
state is known to be the unique ground state of a local, 
gapped graph Hamiltonian. By local we mean involving 
direct interactions among only a fixed number of par- 
ticles; and by gapped we refer to iV-independent finite 
energy difference between the ground and first-excited 
states. These Hamiltonians make adiabatic state cre- 
ation possible: By engineering the interactions so as to 
effectively reproduce the Hamiltonian, the state can be 
prepared simply by first cooling down the system to zero 
temperature and then switching on the interactions so 
that the system is adiabatically driven toward the ground 
state of the graph Hamiltonian. The energy gap imposes 
a threshold to the energy that the environment or erro- 
neous manipulations have to pump into the system in 
order to drive it out of the ground state. This pecu- 
liarity provides the adiabatic approach with an intrinsic 
robustness that the operational approach does not pos- 
sess. In addition, the adiabatic approach is naturally 
better suited for state preparations at large scales. This 
explains the great theoretical effort devoted to finding 
local gapped Hamiltonians with universal resources for 
MBQC as ground states in the DV case, where consider- 
able progress has been achieved [J- On the other hand, 
no such Hamiltonian has been reported for the CV do- 
main. This is what this Letter presents. 

We derive CV graph Hamiltonians Hq that are 
quadratic, gapped, frustration-free (defined below), and 
of range two (coupling only nearest or next-to-nearest 
neighbors). The gap is independent on N and is propor- 
tional to the squared inverse of a squeezing parameter 
s. The ground states are the CV Gaussian graph states 
[2 HI El E], which become universal resources for MBQC 
[2] at criticality s — > oo. With Gaussian states as ground 
states, the fact that Hg is quadratic — and therefore two- 
body — is no surprise. However, this is still in striking 



contrast to the case of qubits, where it is known that no 
state useful for MBQC can be the unique ground state 
of any two-body frustration-free Hamiltonian [10 . 

With Hg at hand, it is now also possible to assess the 
properties of these bosonic systems in the realistic case of 
thermal equilibrium at nonzero temperatures. For exam- 
ple, we show that the correlations across any multiparti- 
tion of the total thermal state are determined exactly by 
those of the boundary subsystem in a thermal state at 
the same temperature. By correlations, we refer to those 
measured by any quantifier invariant under local unitary 
transformations, entanglement being arguably the most 
prominent example thereof. The boundary subsystem 
is composed by the bosons lying at the boundary of the 
multipartition, and is typically much smaller than the to- 
tal system. So, a considerable decrease in computational 
effort is gained for the correlations calculation. In turn, 
this automatically delivers a correlation area law |llj . 
even at criticality. With this, in addition, we prove that 
for any s < oo, these systems typically display thermal 
bound entanglement, even in the thermodynamic limit 
N ^ oo. 

CV graph Hamiltonians. — Let us start by recalling the 
operational definition of graph states for CV quantum 
modes (qumodes) [21 14], which is completely analogous 
to that of the DV case: First, for each j G V, initial- 
ize qumode j in the zero-eigenvalue eigenstate |0)pj of 
the momentum-quadrature operator pj. The computa- 
tional basis is taken as that of the eigenstates \vj)qj 
of the position-quadrature operator g^, with [qj^pk] = 
i^jk,^ j,k G V (we take h = I throughout). The two 
bases are related via the Fourier transform Fji \vj)p. = 
■^ J^duje'-'^^'"^\uj)q^ = Fj\vj)q. . Our starting point 
|0)p^ is thus the uniform superposition of all computa- 
tional states, exactly as in the DV case. Second, for each 
{j, k} € £, apply a maximally-entangling controUed-Z 
gate CZjj. = e^''^^^'' on neighboring qumodes j and k. 
The resulting state is the CV graph state 



\G)=CZ\0)p, 



(1) 



where CZ (without subindices) is a short-hand notation 
fo^ I\{j,k}e£ '^^jk and \v)p = <S)jev I'"j)pj' ^it^ v (with- 
out subindex) standing for the multi- index {vi, ... fjv) G 
M^. In Eq. ^ (and throughout the Letter) we use 
"G" to represent "G(v_£)", unless explicitly indicated. 
Finally, it is also convenient to introduce the nullifiers 



N,, 



Pi 



je^f^ 



(2) 



with Afi all the neighbors of i, whose null-eigenvalue 
eigenstates are the CV graph states: Ni\G) =0, Vi G V 
[2]. The latter is a necessary and sufficient condition to 
univocally specify state ([I]). 

Momentum eigenstate |0)p. can in turn be obtained 
by infinitely squeezing the vacuum coherent state |0)j: 



\0)p. = Muis-yoo Sj{s)\0)j. The action of the uni- 
tary squeezing operator Sj{s) = e*'"(''^'^'^jPj"'"P^'3j)/2 jg ^q 
squeeze the position quadrature by a factor of s and 
to stretch the conjugate momentum quadrature by a 
factor of 1/s: Sj{s)qjSj{s) = qjS and SJ{s)pjSj{s) = 
Pj/s. Thus, states (fTl) can be obtained from the 
vacuo |0) = {S);gv \^h ^^ ^^^ following way: \G) — 
lims^ocCZS{s)\0), being S{s) = {S>iev '^i('^)' "^^^^ * = 
(si, ... sn). In a more general way, finitely-squeezed 
Gaussian graph states are defined as 



\Gs)^U{s)\0) = CZS{s)\Q), 



(3) 



where U{s) = GZS{s) has been introduced. 

Now, consider the eigen-equation of N free, non- 
interacting harmonic oscillators in the ground state. 



ifo|0)^5]^((zf+pn|0)=So|0), 



(4) 



iev 



with Hamiltonian Hq, angular frequencies cui > 0, and 
ground-state energy Eq = X^iev'^i/^- Next apply the 
unitary operator U{s) to both sides of Q from the left: 



f^(^)E ?('?' +pf}UHs)U{s)\0) ^ i?o|G, 



(5) 



iGV 



Invoking the quadrature transformations under squeez- 
ing above we can rewrite the last equation as 
C^^E,6vf(9'A'+pfs2)GZt |G,) = i?o|G,), where 
S'^{s) = 5'(l/s) has been used. The remaining controUed- 
Z gates commute with all position operators but trans- 
form each momentum operator as CZj.pjGZij = pi + qj. 

Thus it is CZpjCZ'^ = Nj. Equation ^ then takes the 

form Y.^ev f (^'Z^' + ^'s')|G.) = ^o|G,), which can 
be renormalized to the more convenient form 



Hg{s)\Gs) = ^\Gs), with 



Ho{s)^J2^{qVs' + N^). 



iGV 



(6) 
(7) 



Equation ([6|) constitutes a new ground-state eigen equa- 
tion, with (ITI) as the new Hamiltonian, Q as the new 
ground state, and —§- as the new ground-state energy. 

Hamiltonian ([7]) is in turn the desired graph Hamilto- 
nian. The natural appearance of the nuUifier operators 
in it is remarkable. As anticipated, the two-body char- 
acter of Hg{s) is encapsulated in its quadratic form. In 
view of Eq. (l2|, it is clear that direct couplings only 
take place between nearest, or next-to-nearest, neigh- 
bors. In addition, each and all of the terms in (It]) com- 
mute, which implies that the Hamiltonian is frustration- 
free (meaning that the ground state minimizes the en- 
ergy of each local term in the sum). Furthermore, the 
gap between the ground and first-excited states is that 
of Hq (wi^jn = minr-jgyi oji) consistently renormalized: 



ujimin/s^ ■ At infinite squeezing the gap vanishes and the 
system is then called critical. Interestingly enough how- 
ever, its thermal states always satisfy a correlation area 
law, as we show below. In addition, at criticality the 
graph Hamiltonian ([7]) acquires an even simpler form: 

Hg = lim.^oo Hg{s) = E,ev ^^f ■ 

The derivation of ^ not only delivers the desired 
Hamiltonian but also comes with the interesting byprod- 
uct of readily giving the symplectic transformation that 
takes the Hamiltonian to its normal mode decomposition. 
In this case, U{s) is the unitary representation of such 
transformation and maps our Hamiltonian to that of a 
collection of noninteracting harmonic oscillators: 

U\s)Hg{s)U{s) = S\s)CZ^Hg{s)CZS{s) = Ho/s'^,{8) 

with Ho given in Q. What is more, the same uni- 
tary transformation delivers also the Gaussian nulli- 
fiers. That is, N commuting operators with the Gaus- 
sian graph state ^ as their (unique) mutual eigenstate 
of null eigenvalue. To see this, instead of eigen equation 
Q, start by aj\0) = 0, with Oj = (qj + ipj)/^/2 the 
annihilation operator of the jth qumode, and with the 
same reasoning as above arrive at the Gaussian nuUifiers 
Nj{s) = -iqj/s^ + Nj, where Nj is the jth nuUifier ^. 
Thermal Gaussian graph states. — Once we have ob- 
tained Hamiltonian (l7| we can now consider thermal 
Gaussian graph states, which are defined in the usual 
way: 



-Hg/T 



PG,T 



Tr[e 



-Hc/n 



(9) 



where T is the temperature of the system's thermal bath 
(Boltzmann's constant is set as unit). 

We are now in a position to study the correlations 
C{pg,t) across any multipartition of thermal state Q, 
with C any arbitrary correlation quantifier invariant un- 
der local unitary transformations. For its evaluation 
we first decompose the symplectic unitary as U{s) = 
CZS{s) = CZxSy{s) (g) CZ^Syis) = Uy{s) (g) Uy{s). 
A" C £^ is the set of boundary-crossing edges [12 , those 
that connect vertices belonging to different subpartitions. 
The latter vertices in turn compose the set 3^ C V of 
boundary vertices jl2| . The two sets together consti- 
tute the boundary subgraph G(^y,x)i a-nd the rest G/tj^n, 
with y = V/y and X = £/X, is called the non- 
boundary subgraph. With this, we notice that pg,t = 
Uyis)pG,y,;^^,T «) pO(y),Tt/^(s). In the last, PG(y.;t),T is 
a thermal state of the boundary subsystem, defined as 
in Eq. ^ but with respect to the boundary subgraph 
Hamiltonian Hc^y,,^ = 'E.ey fi^f/^'^ + ^yz)- Here, 
Nyi = Pi — J2i£j\f i Qj is the ith nuUifier corresponding 
to the boundary subgraph - the same as in Eq. (pi) but 
with the sum running over the set Myi = Afi oy oi the 
neighbors of i in 3^ -. In turn, Po,-, ,t is a thermal state of 




FIG. 1: High-temperature, high-squeezing region of the dis- 
tillability phase diagram of thermal Gaussian graph states 
associated to arbitrarily large linear clusters. The solid 
curve represents the inverse critical temperature Tg^^ = 
2 arccoth(\/l + s^ + s^) (in arbitrary units) at which the 
negativity between any two contiguous blocks vanishes as a 
function of the inverse squeezing, and the dashed one the anal- 
ogous Ts"^ = 2 arccoth[(\/l + 4s4 - 2^2^28^ + s"*)"^] for 
the bipartitions of any qumode versus the rest. Between the 
two curves every qumode is entangled with all the rest but no 
entanglement can be locally extracted. 



the nonboundary subsystem with respect to the decou- 
pled harmonic Hamiltonian Hoy = J^iey '^i'li + Pi)- 
Notice that HGi^y ^^ and Hq commute. 

Now, by definition, Uy{s) is a local unitary opera- 
tion with respect to the considered multipartition, so 
we can disregard it as for what the correlation eval- 
uation concerns. Once Uy{s) is omitted, the bound- 
ary and nonboundary subsystems are left in the prod- 
uct pG^y x),T ® Po,-yT, so all the correlations across the 



multipartition are concentrated in its boundary: 



C{pg,t) = C{pc 



^:T 



(10) 



Thermal states of one-dimensional bosonic chains gov- 
erned by some specific families of finite-ranged quadratic 
Hamiltonians are known to satisfy an entanglement area 
law for some particular entanglement quantifiers [llj . 
An area law is said to be satisfied when the correla- 
tions across a multiparition scale at most with the size 



of its boundary. Equivalence ( 10 1 gives us the fully gen- 



eral statement for thermal Gaussian graph states: they 
obey an area law for any geometry and any local-unitary 
invariant correlation. Not only that, it gives us much 
more refined information for it reduces the correlation- 
evaluation problem of arbitrarily sized specimens (for ex- 
ample macroscopic ones) to one on the boundary subsys- 
tem, regardless of how well an area law is satisfied. 

In Fig. [TJfor instance we have plotted the inverse crit- 
ical temperatures at which the logarithmic negativities 
[13] of state (l9| vanish as a function of the inverse squeez- 
ing, for the exemplary case of one-dimensional graphs 



and constant Hamiltonian couplings uji = tu. The solid 
curve corresponds to any block-to-block bipartition, and 
the dashed one to any one-mode-versus-the-rest biparti- 
tion [see inset of Fig. fT|. By equivalence (10 1, the former 
is given simply by the critical temperature T2c of a two- 
qumode thermal cluster state, and the latter by that of 
the bipartition of one qumode versus the other two in 
a three-qumode thermal linear-cluster state, T^^.. Both 
temperatures can be calculated analytically. For temper- 
atures higher than T2c not even two-mode entanglement 
can be (locally) distilled between any two qumodes in the 
chain, for if this were possible then there would be two 
contiguous blocks with positive negativity [HI [15] . On 
the other hand, for temperatures lower than (at least) Ta^ 
every boson in the chain is entangled with all the rest. 
For T2c < T < Tsc thus, thermalization naturally drives 
Gaussian 1 D graph states of any size to bound entan- 
glement. This extends of course to higher-dimensional 
clusters and nonequal couplings, as in the DV case |14) . 

Phase-space diffusions on CV graph states. — In the 
limit of infinite squeezing, thermal states Q become 
equivalent to the result of independent Gaussian diffu- 
sion along the q direction on pure states ([I]). This in- 
teresting connection between (collective) thermalization 
and (independent) dephasing is the CV version of the 
one observed in qubit graph states [UJ [TC] . Additionally, 
also for s — )• oo, the evolution of correlations under noise 
processes described by arbitrary phase-space-shift maps 
can be monitored in terms of the boundary subsystem 
by translating the whole machinery for the study of DV 
graph-state entanglement under Pauli maps developed in 
Ref. jl2j. This characterization is relevant to the devel- 
opment of CV quantum error correction schemes, but it 
will be touched upon elsewhere. 

Discussions. — As known, states (ll]) are as nonphysi- 
cal idealizations as the free-particle states |0)p. However, 
universality for MBQC has only been proven for infinitely 
squeezed states ^. What is more, recent results j9j show 
that large-scale MBQC with "imperfect" physical states 
([3]) faces fundamental limitations that may only be cir- 
cumvented with the full machinery of quantum error cor- 
rection and fault tolerance [6 . The latter is yet to be de- 
veloped for CV systems, but will presumably demand a 
very large overhead in resources. This in turn highlights 
the importance of local, short-ranged, gapped CV graph 
Hamiltonians, for it is precisely the large-scale scenario 
where the adiabatic approach is specially well suited. On 
the other hand, states pi do yield universal resources for 
small-sized computations. 

From a more applied viewpoint, our findings open 
new realistic venues for the physical realization of CV 
quantum computing beyond the standard optical ap- 
proaches. In fact, the basic constituents of Hamiltonian 
(It]) — namely, the couplings qi ® qj and pi ig) qj — have 
already been demonstrated in technologically mature ex- 



perimental platforms, as Coulomb crystals [TTj and opto- 
mechanical resonators [18] . In addition, this type of cou- 
plings have also been envisioned in coupled-microcavities 
arrays or superconducting waveguides |19| . All these con- 
stitute examples of versatile and promising architectures 
where controlled geometrical arrangements of the desired 
interactions seem feasible in a very near future. 

Note added. — After conclusion of this work, another 
paper addressing CV graph Hamiltonians appeared [20] . 
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